Acoustic decoherence of flux qubits 
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Decoherence of a flux qubit due to inelastic scattering of thermal phonons by the qubit is studied. 
The computed decoherence rates contain no unknown constants and are expressed entirely in terms 
of measurable parameters of the qubit. The answer depends strongly on the size of the qubit as 
compared to the wavelength of phonons of frequency / = A/(2irh), with A being the tunnel splitting 
of the qubit. Thermal phonons set the upper limit on the operating temperature of a small qubit 
at around 10 — 20K. For large qubits acoustic decoherence due to one phonon processes should be 
taken into account at any temperature when a high quality factor is desired. 
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I. INTRODUCTION. 



Flux qubits have received much attention lately due 
to low decoherence exhibited at a macrosopic scale. In 
essence such a qubit consists of a small superconducting 
loop interrupted by one or more Josephson junctions. 
When half a flux quantum is applied to the loop, the 
system forms a symmetric double-well potential. The 
two classical states associated with the minima of this 
potential correspond to clockwise and counterclockwise 
circulating currents. The fascinating property of such a 
system is that it exhibitsquantum tunneling between the 
classical minima [H, HI HI 0]' thus, providing an example of 
a quantum superposition of macroscopic current states. 
The ground state and the first excited state of the sys- 
tem are symmetric and antisymmetric superpositions of 
clockwise and counterclockwise current states, separated 
by the tunnel splitting A. When the system is prepared, 
e.g., in a clockwise current state, the probability to find it 
with clockwise or counterclockwise current oscillates co- 
herently at a frequency / = A/(2ivh). Experiments have 
demonstrated that the quality factor of the flux qubit 
can be as high as 5000 Q, thus, pushing it to the front 
line of promising candidates for quantum computation. 
Recently, this promise has been further amplified by suc- 
cessful experiments with coupled flux qubits @, 0, @] • 

Most of the work on decoherence of flux qubits concen- 
trates on the effects of non-zero impedanc e 1131 . electro- 
magnetic 0, H OH EH and 1// noise pj, ^7 Decoher- 
ence due to coupling of a flux qubit to the acoustics waves 
has been studied to a lesser degree @. Meantime, such 
a decoherence is generic for a flux qubit. This can be 
seen from the following argument. Tunnelling between 
clockwise and counterclockwise circulation is accompa- 
nied by the reversal of the angular momentum associated 
with the circular current. In experiments on flux qubits, 
this angular momentum can be as large as 10 10 ?i, so that 
its non-conservation would be quite dramatic. A similar 
problem exists for any LC circut. To conserve the angular 
momentum, a freely suspended ring containing an induc- 
tor and a capacitor would have to co- wiggle mechanically 
with the oscillating current. When such a ring is attached 



to a solid matrix it should produce torsional oscillations 
in the matrix. This immediately suggests that unless the 
flux qubit is in the ground state, it should have a finite 
probability to radiate a phonon. The latter can be viewed 
as a consequence of the fact that the electric current is 
defined with respect to the crystal lattice of ions so that 
quantum states of the current are inevitably entangled 
with the lattice states. 

It was shown earlier Q that the effective coupling of 
the oscillating current with the crystal lattice has an uni- 
versal form that does not contain any fitting parameters. 
At very low temperatures the typical decoherence rate 
due to this coupling ranges from Is -1 for a small qubit 
to 10 6 s _1 for a large qubit. The low temperature deco- 
herence rate is due to a spontaneous decay of an excited 
current state to the lower energy state, accompanied by 
the emission of a phonon. Such processes are known as 
one phonon processes. They dominate relaxation and de- 
coherence due to phonons in the low temperature limit. 
Since operating qubits at elevated temperatures is de- 
sirable one must consider other phonon processes which 
may be dominant at these temperatures. Such processes 
consist of inelastic scattering of thermal phonons by the 
qubit. They are known as two phonon processes. In this 
paper, we focus on the decoherence via a two phonon 
channel. We find that the significance of this channel 
strongly depends on the size of the qubit, R, as com- 
pared to the wavelength, A = ft//, of a phonon of fre- 
quency / = A/(2ttH) (with vt being the speed of the 
transverse sound). It turns out that for a small qubit 
the two phonon processes can be the main factor limit- 
ing the qubit operating temperature. Another side of the 
coin, however, is that the corresponding two phonon rate 
is proportional to the square of the energy bias between 
the two states with opposite directions of the circulating 
current. This, in principle, allows one to switch the two 
phonon decoherence on and off by changing the external 
magnetic flux that controls the bias. For a large qubit 
the one phonon processes seem to always win over the 
two-phonon processes. 

The structure of the paper is as follows. Interaction 
between superconducting current and elastic twists is dis- 
cussed in Section II. Two phonon processes are studied in 
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Section III. Subsection III-A deals with inelastic scatter- 
ing of thermal phonons by small qubits. Matrix elements 
and decoherence rates for large qubits are studied in Sec- 
tion III-B. Section IV contains comparison of one-phonon 
and two-phonon rates, as well as numerical estimates of 
the acoustic decoherence rate, and general discussion of 
the results. 



II. COUPLING OF SQUID STATES TO 
ACOUSTIC PHONONS 



We consider a loop of radius R carrying a current J 
oscillating with frequency /o ~ 10 9 — 10 10 s _1 between the 
clockwise and counterclockwise direction of motion. The 
angular momentum associated with the circular current 
changes from L to — L every half a period. As a conse- 
quence, the local environment co-wiggles with the same 
frequency and produces elastic distortions of the lattice in 
accordance with angular momentum conservation. These 
distortions, u(r,i), being pure twists, satisfy 



where 



0. 



(1) 



The kinetic energy associated with the superconduct- 
ing current can be written as 



KE 



: (Viat + u) 2 , 



(2) 



where v e is the velocity field of electrons in the labora- 
tory coordinate frame, V£ Qt is same velocity in the lat- 
tice frame, n e is the concentration of electrons, and m e 
is electron mass. Interaction of the current with lattice 
distortions can be formulated by noting that the current 
is defined in the frame of the moving ions, 

j = en e w Lat = en e {w e - u), (3) 

where e is electron charge. From Eq. ([2|) one obtains Q : 



'Hint = — / crrj • ii. 



(4) 



This Hamiltonian satisfies all the symmetries and is free 
of any unknown interaction constants. 

We denote the eigenstates of the angular momentum, 
L z — L and L z — —L, as | f) and | J.), respectively. 
The system under consideration can be approximately 
modeled as a particle of spin L in a biased double well 
potential described by the Hamiltonian H.q whose lowest 
eigenstates are 



10) 



V2 



(C_|T) + C + | j» 



and 



Here 



_(C + |T>-C_|i». 



C ± = Vlie/A, 



(5) 



(6) 



(7) 



A = 



A 2 , 



(8) 



is the energy splitting of |0) and |1), e is the bias, and 
A is the energy splitting at £ = 0. 

The total Hamiltonian of the system is 



H = Ho + H p h + Hi 



(9) 



The first two terms stand for the interaction-free qubit 
and phonon Hamiltonian, respectively, and the last term 
is given by Eq. Decoherence of a general superpo- 
sition state |\&) = ci|0) + C2 1 1) will take place through 
relaxation of |1) to |0). The corresponding two phonon 
transition is determined by the matrix element of Hint 
that will be computed in the next Section. To shorten 
formulas, unless stated otherwise, we will use the units 
in which fcs = fi. = 1. 



III. TWO-PHONON PROCESSES 

A. Small qubit. 

In the case of a small flux qubit of size R <C A, one can 
treat the local environment at the position of the qubit as 
a rigid matrix making a tiny local rotation with angular 
velocity (13 ] 



= -V x ii. 
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(10) 



Equivalently, the velocity vector u can be expressed in 
terms of the angular velocity as ii = SI X r, where r is 
the position vector with its origin at the center of the 
qubit. The interaction Hamiltonian of Eq. ^ can now 
be simplified as 



Hint — L • n, 



in which 



L = ^ 



d r r x j 



(11) 



(12) 



is the angular momentum of the circulating current. 
Hamiltonian of Eq. (fTTj) is a consequence of the con- 
servation of angular momentum. 

We are interested in the process of relaxation through a 
two phonon channel in which one phonon is absorbed and 
another is emitted with wavevectors k and q respectively. 
The matrix element for this process is 



M 



E 
E 



{^ f \n in t\^) mnint^i 

(V f\H int \lps) (^\Hint\% 



Ex 



E : 



(13) 



where l^i) (|^/)) stand for the initial (final) state of the 
system, and w q are phonon frequencies, and E\ is 
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the energy of the first excited state |1). The summation 
over £ labels the energy states of Ho, and the interme- 
diate phonon states |nk — l,n q ) in the first term and 
|?Jk,«q + 1) in the second term. With the help of Eq. 
(|10P and conventional quantization of the phonon field, 



u 



— y 



(«kA + al kjA ) - ( 14 ) 



where p is the mass density of the solid matrix and ekA 
is the polarization vector of the phonon, one obtains for 
the phonon part of the matrix element 

M p h = (n q + 1| Q z |n q ) (rc.k - 1| 0, z |nk) 
= (ra k - 1| fi, |n k ) (rc. q + 1| Vl z |n q ) 

= ^f/I k x e x } z [q x e cr ] 2 y / w k A^ qC r("- q + l)"k- 

(15) 

The action of the angular momentum operator on |0) and 
|1) produces the following states: 



L z \l) = 



LAo 

A 
LAp 

A 



|0> 



A 

e 



|0) 



A 



ID 







(16) 



Inserting these results into Eq. (fT5|) one can immediately 
see that \tp^} must be either |0) or |1). Thus, we only need 
to consider 



-<0|L Z |0) = (1|L,|1) = ^ 



\L Z \1) 



(1\L, 



LA 



(17) 



Taking into account the conservation of energy, w q — cjk 
A, the matrix element (IT31) reduces to 



M = M, 



2L 2 eAo 



1 



ph ' 



+ A) 



(18) 



The relaxation rate can be obtained using the Fermi 
golden rule in the second order, given by 



r? 



E 

A.CT 



dkdq 



V 2 \M\ 2 2TrS(ujq — cok— A). (19) 



Substitution of Eq. (18]) into Eq. (HD yields 



r4_2 A 2 

16p 2 A 2 (2tt) 5 



J dkdq[(k x e A ) z (q x e a ) z f 



Wk(wk + A) 2 



S(0Jq - CJk - A). 



(20) 



Replacing the integration variable <ik with oj^duj^dil^ / vf 
and the vector k with kwk/fi and k = k/fc, one obtains 

r 2 = 



L 4 £ 2 A 2 B 2 



16p 2 A 2 (2tt) £ 



duui 3 (uj + A) 3 



(e^-l)(l-e-("+ A )/ T )' 



(21) 



where 



B = E/^xeA)^. 



(22) 



The sum in Eq. (|22p runs over the two transverse polar- 
izations. Making the change of variable x = lu/T yields 



where 



L^AlB^ 
T * = 16(^)W ff( ' ^ 

dxx 3 (x + A/T) 3 



fl(A/T) = / 
Jo 



(e x -l)(l-e- x - A / T )' 



(23) 



(24) 



For completeness, we insert kg and /i into our result 
which in the limit fcgTc A reduces to 



7T L 4 £ 2 A 2 Afc4T 4 



(25) 



270 p 2 t; 10 ?i 7 
In the opposite limit fc^T ^> A one obtains for the rate 
_ _ 15tt L 4 e 2 A 2 fc 7 T 7 



36 p 2 A 2 vj°h 7 



B. Large qubit. 



(26) 



For a large qubit, i? 3> A, one must integrate the in- 
teraction of the superconducting current with phonons 
along the entire loop as in Eq. (0J , 



Hint 



Wk 



— E\/^Bjk-eOKA-atkA), ( 2? ) 



where jk = / dr 3 je lkr is the Fourier transform of j. The 
solution to jk in the thin-ring approximation is 



j k = -i2TrRJ 1 (k ± ,R)Jn k , 



(28) 



where 1 k is a unit vector in the plane of the ring, 
is the z-component of k, and J\ , R) stands for 
the Bessel function of the first order. The interaction 
Hamiltonian can now be written as 



L • Sl eff , 



(29) 



where 

Meff = 



R \ 



Ji(k±,R)(a,kt - afkt)e z - (30) 



The scalar product nk • e A = S t \ as nk is in the plane of 
the ring. To obtain the two phonon matrix element one 
needs only to replace (k x e A ) z (q x e a ) z in Eq. (TT5"|) with 



Jl(k±R)Ji(q±R). 



(31) 
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Then, the rate for a large qubit becomes 

TtL 4 £ 2 A 2 R 4 



r 2 



j-l-j- / dw/(u,A > r)e(w J ii,«t,A) ) 



where 



/(w,A,T) 



w(w + A) 



( e w/T _ _ e -(w+A)/T): 



(32) 
(33) 



and 



®(gj,R, ft, A) = / d6>isin6>iJ 2 
Jo 

/•7T 

x / d8 2 sin 9 2 J 2 
Jo 



luR 
v t 



■ smtq 



(34) 



Due to its large argument, Ji([u) + A]R/vt sin0 2 ) can be 
approximated by its asymptotic form. The integral over 
9 2 in Eq. (13~4")) then yields approximately «t/(w + A)i?. 
The integral over 6\ equals 



1 fojR\ 2 /3 5„ /wi?\ 2N 



3 V t; t / 2 2 \ v t 



(35) 



where i-F^ is the generalized hypergeometric function. In- 
serting the results just obtained into Eq. (|32|) and making 
the change of variable x — lu/T, one obtains 

^ irL 4 e 2 A 2 R 5 T 4 p /T x 3 (e x - l)- 1 i J F 2 (-/3 2 



24p 2 A 2 a 4 vJ J 



(i.g-x-A/T) 



(36) 

where /3 = RT/vt- In the limit T» A, the integral in 
Eq. (|36|) simplifies to 



00 x 3 e x 1 F 2 (-f3 2 x 2 ) 
(e* - l) 2 • 



(37) 



or, after integrating by parts 



*^ A ':y»4r^. (38 ) 



/? 2 7 



e x - 1 



For f3 > 1, Eq. (H§ reduces to 3/3 -2 . The rate then 
becomes 



„ 7rX 4 £ 2 A^ 3 fc|T 2 
2 ~ 8p 2 A 2 aVft 2 ' 



(39) 



where we again inserted /cb and /i into our result. In the 
limit feT C A, one obtains 



r 2 s* 



nL 4 e 2 A 2 R b k%T 4 
2Ap 2 A 2 a 4 v 7 t h 4 



K 



Rk B T 
hv t 



where 



(40) 



(41) 



IV. DISCUSSION 



We are interested in the transition temperature be- 
tween the rate due to a direct process and the Raman 
process. The former was calculated in Ref. 0. It was 
shown that 



L 2 A 5 , / A 
■ coth 



12npvfh 4 """" \2k B T 
for a small qubit and 

L 2 A 2 , / A 



Ti = 



4irR 3 pv 2 h 



coth 



2k B T 



(42) 



(43) 



for a large qubit. 

In the limit k B T < A, equations ((25j) and (|42j) yield 
for a small qubit 



£2 _ 12ir 2 L 2 e 2 A 2 k 4 B T 4 
fT " 



270pA 4 vfh 3 



(44) 



For, e.g., L ~ 10 2 , p ~5g/cm 3 , j) t ~5x 10 3 m/s, and / ~ 
5 x 10 9 s- x , the ratio in Eq. (05) yields 10" 8 T 4 . Clearly, 
for k B T <C A (T ~ 10~ 2 K), the two phonon process is 
utterly insignificant. In the opposite limit k B T 3> A, 
equations (p6|) and (|42|) give 



(45) 



£2 _ 907r 2 L 2 £ 2 A 2 fc|r 6 
Ti ~~ pA & v b t h 3 



For the same parameters as considered above, Eq. 
yields 10~ 2 T 6 . The two phonon process thus dominates 
over a one phonon process above T ~ 2K. At T ~ 20K 
the decoherence rate due to two phonon processes is of 
order 10 6 s _1 and, therefore, its contribution to the qual- 
ity factor of the qubit cannot be ignored. 

For a large qubit in the small temperature limit one 
obtains with the help of equations (|40|) and (|43| 



= L 2 e 2 A 2 k%T 4 ( Rk B T \ 
T 1 6ir 2 pA 4 vfh 3 \ hv t )' y ' 

For a large qubit parameters: L = 10 10 , R ~ 10 _4 m and 
T ~ KT 2 K, Eq. (gH) yields 1(T 5 . As in the case of a 
small qubit, the two phonon rate is negligible compared 
to the one phonon rate. In the large temperature limit 
equations (f3"9"| and (|13"|) give 



L 2 e 2 A 2 fc s T 



Ti 4 7 r 2 i? 2 j oA 3 i; 3 ft' 



(47) 



which yields 10 -6 T. Evidently, for a large qubit, one 
phonon processes are the dominant source of decoherence 
at any temperature. For the numbers used above the one 
phonon rate is of order 10 6 s _1 . 

A few important observations for flux qubits follow 
from the above results. Firstly, acoustic decoherence 
should definitely be taken into account when the qual- 
ity factor as large as 10 4 is desired. Secondly, for small 
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biased qubits the inelastic scattering of thermal phonons V. ACKNOWLEDGEMENTS 

by the qubit can become a dominant mechanism of de- 
coherence above 20K. Finally, the proportionality of this 
mechanism to the bias provides an additional way to con- 
trol the flux qubit. The interesting feature of the above 
results is that they do not contain any unknown constants 

- the decoherence rate is expressed entirely in terms of This work has been supported by the Department of 
measurable parameters of the qubit. Energy through Grant No. DE-FG02-93ER45487. 
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